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We use an exact renormalization-group transformation to study the Ising model on a complex 
network composed of tightly-knit communities nested hierarchically with the fractal scaling recently 
discovered in a variety of real- world networks. Varying the ratio Kj J of of inter- to intra-community 
couplings, we obtain an unusual phase diagram: at high temperatures or small K/ J we have a 
disordered phase with a Griffiths singularity in the free energy, due to the presence of rare large 
clusters, which we analyze through the Yang-Lee zeros in the complex magnetic field plane. As 
the temperature is lowered, true long-range order is not seen, but there is a transition to algebraic 
order, where pair correlations have power-law decay with distance, reminiscent of the XY model. 
The transition is infinite-order at small K/J, and becomes second-order above a threshold value 
(K/J) m . The existence of such slowly decaying correlations is unexpected in a fat-tailed scale-free 
network, where correlations longer than nearest-neighbor are typically suppressed. 
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I. INTRODUCTION 

Many real-world networks have modular structure [J 
& IE 0] : their nodes are organized into tightly- knit com- 
munities where the node-node connections are dense, 
with sparser connections in-between communities. This 
structure is often hierarchically nested, with groups of 
communities themselves organized into higher level mod- 
ules. Given the relevance of modularity to features like 
functional units in metabolic networks Q, it is not sur- 
prising that community structure has become one of the 
most intensely studied aspects of complex networks. Re- 
cently, Song, Havlin and Makse [E [(| discovered that 
certain modular networks also possess another remark- 
able characteristic: fractal scaling, where the hierarchy of 
modules shows a self-similar nesting at all length scales. 
Examples of such fractal networks include the WWW, 
the actor collaboration network, protein interaction net- 
works in E. coli, yeast, and humans, the metabolic path- 
ways in a wide variety of organisms [H, and genetic reg- 
ulatory networks in S. cerevisiae and E. coli [7J. 

Despite the widespread occurrence of fractal topolo- 
gies, little is yet known about the nature of cooperative 
behavior on these networks, or even more generally on 
how modular structure affects collective ordering or cor- 
relations among interacting objects. In particular the 
Ising model has been investigated extensively on non- 
fractal scale-free networks @, S, HE GH, H, HE M, HE EH , 
but only recently has a form of community structure 
been included: an Ising ferromagnet was studied on two 
weakly coupled Barabasi- Albert scale-free networks with 
a varying density of inter- network links [ItJ , finding sta- 
ble parallel and antiparallel orderings of the two com- 
munities at low temperatures. It would be interesting 
to examine a system with a large number of interacting 
communities, capturing more fully the complex modu- 
lar organization of real-world examples. In this paper, 
we introduce a hierarchical lattice [IE EE HE] network 



model exhibiting a nested modular structure with frac- 
tal scaling. Hierarchical lattices (part of a broader class 
of deterministically constructed networks [2l|, HE HE HH 
HE HE HE HI Hi HE]) have been the focus of increas- 
ing attention recently [3l|, HE HE H3], since they can 
be tailored to exhibit various features — including scale- 
free degree distributions, small-world behavior, and frac- 
tal structure — for which exact analytical expressions can 
be derived. The explicit results from such determinis- 
tic models can serve as a testing ground for approximate 
phenomenological approaches, and a starting point for 
extensions incorporating additional realistic features like 
randomness [34j |. 

Here we exploit another advantage of such lattices: the 
ferromagnetic Ising model can be solved through an ex- 
act renormalization-group (RG) transformation. Varying 
the strength of interactions between communities, we find 
an unusual combination of thermodynamic properties. 
At high temperatures or weak inter-community coupling 
the system is disordered, but the free energy as a func- 
tion of magnetic field H is nonanalytic at H = 0. This 
is due to the presence of rare large communities, similar 
to the Griffiths singularity in bond-diluted ferromagncts 
below the percolation threshold [35| : there the system is 
partitioned into disjoint clusters of connected sites, and 
the small probability of arbitrarily large clusters leads 
to an analogous nonanalyticity in the free energy above 
T c . As we lower the temperature in our network, true 
long-range order is never achieved at T > 0, even for 
the strongest inter-community coupling. Surprisingly, we 
find instead a low-tcmperature phase with algebraic or- 
der, just as in the XY model: the magnetization is zero, 
but there is power-law decay of pair correlations with dis- 
tance, and the thermodynamic functions throughout the 
entire phase behave as if at a critical point. 

The organization of the paper is as follows: in Sec. II 
we describe the network's construction (Sec. II. A) and 
summarize its topological properties (Sec. II. B), includ- 
ing its community structure and fractal scaling charac- 
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FIG. 1: (a) Construction of the hierarchical lattice, (b) The 
lattice after t — 4 steps in the construction. Communities at 
the n = 0, 1, 2 levels of the hierarchy are shown with white, 
dark gray, and light gray backgrounds respectively. 



teristics. Sec. Ill examines thermodynamic properties of 
the Ising model on the network, derived from an exact 
renormalization-group approach (Sec. III. A). We discuss 
the phase diagram and critical behavior in Sec. III.B, and 
then focus on two particularly interesting aspects of the 
results: the presence of Griffiths singularities in the free 
energy (Sec. III.C), and the nature of long-range pair cor- 
relations in the low-temperature phase (Sec. III.D). We 
present our conclusions in Sec. IV, and note that the be- 
havior described here is characteristic of a broader class 
of hierarchical lattice complex networks — a fact that will 
be explored in future studies. 



II. NETWORK PROPERTIES 

A. Construction procedure 

Our lattice has two types of bonds, depicted as solid 
and dashed lines respectively. At each construction step, 
every solid bond is replaced by the connected cluster on 
the right of Fig. [2(a), and this procedure is iterated t 
times. The initial t = lattice is two sites connected by 
a single solid bond. Fig.QJb) shows the network at t = 4. 
All results quoted below are for the infinite lattice limit, 
t — ► oo. 



B. Topological characteristics 



Size of network: The total number of sites N = | (4* + 
2), the total number of bonds Nj, — |(4* +1 — 1), and the 
diameter of the network (the maximum possible shortest- 
path distance between any two sites) is D = 2 t+1 — 1. 

Degree distribution: The probability P(k) of finding 
a node with degree k is zero except for k = 2 m + 1 for 
some integer m > 1, where P(k) = 3-4~ m . The scale-free 



exponent 7 is calculated from the cumulative distribution 
P cum (k) = £^ =fc P(fc') ~ fc 1 -^. For large k we find 

Pcum(k) f» fc -2 , SO 7 = 3. 

Community structure: We can define several levels of 
hierarchical modular organization, labeled by integer n: 
at the lowest level (n = 0) we have clusters of solid bonds 
(shown with white background in Fig. Hfb)), with the 
dashed bonds acting as inter-community links; at the 
next level (71 = 1) we can group together those com- 
munities which correspond to a single solid-bond cluster 
at the (t — l)th construction step (dark gray background 
in the figure); the n — 2 level communities arc outlined 
in light gray. In general, for a lattice after t construc- 
tion steps, a community at the nth level of the hierarchy 
evolved from a single solid-bond cluster at step t — n. 

Fractal scaling: Adapting the analysis in Refs. 0, 
we can characterize the fractal topology of the network 
through two exponents ds, dk, defined as follows. At 
the nth level, all communities have the same diameter, 
£b = 2 n ^ 2 — 2 (for n < t — 1). Thus at each level the 
communities form a "box covering" of the entire net- 
work with boxes of the same £3- The scaling of the 
number of boxes Nb{£b) required to cover the network 
for a given £b defines the fractal dimension ds, namely 



In our case we have Nb(£b)/N 



N b (£b)/N ~ tg d * 
4-n-i w 4£~ 2 for large n, yielding ds = 2. Similarly 
the degree exponent dk of the boxes is defined through 
&B(^s)Ahub ~ ^B dk , where fcs(^s) is the number of out- 
going links from the box as a whole, and fchub the degree 
of the most connected node inside the box. For boxes 
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with large fchub we get a scaling fcs(£s)//chub ~ 2£ 
giving dk = 1. As with all the real- world fractal networks 
examined in Refs. [1, @, the scale-invariance of the prob- 
ability distribution is related to the fractal scaling of the 
network through the exponent relation 7 = 1 + ds/dk, 
which is satisfied for ds = 2, dk = 1, and 7 = 3. 

Modularity: The strength of community structure — 
the extent to which nodes inside communities are 
more tightly knit than an equivalent random network 
model — is quantified through the modularity [3] Q = 
J2 S [L/Nb — (d s /2Nb) 2 ] , where the sum runs over all 
communities, and l s , d s are the total number of bonds 
and total sum of node degrees for the sth community. In 
our case each level n in the hierarchy describes a different 
partition of the network into communities, and we find 
the corresponding modularity Q — 1 — 4 - " -1 . Thus Q 
increases from 3/4 at n = to the maximum possible 
value 1 as n — * 00, showing that the modular structure 
becomes ever more pronounced as we go to higher levels. 

Distribution of shortest-paths: We define Ni as the to- 
tal number of site pairs (i,j) whose shortest-path dis- 
tance along the lattice 1^ — I. The distance I can take on 
values between 1 and D. Ni has a non-trivial dependence 
on £, but satisfies the scaling form Nt — 2 3t ft(l/ D), 
where the function f t (£/D) approaches for £ close to 
1 or D, and f t (£/D) - O(10" 1 ) for 1< £ < D. The av- 
erage shortest-path length £ — wjfzrn Yli=i ~ D ~ 
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N 1 / 2 , so the network is not small- world. 

III. ISING MODEL ON THE NETWORK 
A. Renormalization-group transformation 

Let us now turn to the Hamiltonian for our system, 

-pH = jJ2 s i s i + K J2 SiS i + H J2 s " (1) 

(ij)s (ij)d i 

where s, = ±1, J, K > 0, and (ij) s , (ij)d denote sums 
over nearest-neighbor pairs on the solid and dashed bonds 
respectively. The ratio of inter- to intra-community cou- 
pling is parametrized by Kj J. The RG transformation is 
the reverse of the construction step: the two center sites 
in every cluster like the one on the right of Fig. HJa) 
are decimated, giving an effective interaction between 
the two remaining sites. The renormalized Hamiltonian 
—(3TC has the same form as Eq. (fTJ), but with interac- 
tion constants J',K',H'. Two additional terms also ap- 
pear: a magnetic field counted along the solid bonds, 
H' B Y2{ij) ( s i + s j')i an d an additive constant per solid 
bond G' . The renormalized interaction constants are 
given by: 

J' = (1/4) In (R^/Rf) , K'=K, H' = H, 
H' B = (1/4) In {Rx/R 2 ) , G' = (1/4) In (R^Rj) , 

(2) 

where: 

Ri = w~ 4 xy~ 2 + 2x~ 1 z 4 + w 4 xy 2 z s , 

R 2 = w~ 4 xy 2 + 2x- 1 z- 4 + w 4 xy- 2 z- 8 , (3) 

i?3 = x w + x w + xy z + xy z , 

and w = e J , x = e , y = e H , z = c Hb . Under renormal- 
ization a nonzero site magnetic field H induces a bond 
magnetic field H B , while H' = H since the site field at 
the edge sites in each cluster is unaffected by the deci- 
mation of the center sites (36j . This transformation is 
exact, preserving the partition function Z' = Z, and 
we iterate it to obtain the RG flows, yielding the global 
phase diagram of the system. Thermodynamic densities, 
corresponding to averages of terms in the Hamiltonian, 
transform under RG according to a conjugate recursion 
relation (36[. Iterating this along the flow trajectories 
until a fixed point is reached, we can directly calculate 
the magnetization M — jj^Z i {si) 1 internal energy per 
site U = jj (H) , and their derivatives x = Tfff > C = ■ 
We are also interested in the pair correlation Gy = 
(siSj) — (si)(sj) for arbitrary sites i, j in the network. 
Since our lattice is highly inhomogencous, Gij is not 
a simple function of However, following the anal- 
ysis in Ref. (37j . we can define an average correlation 
= j):ii-=£} Gij-, where the sum is over all 
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FIG. 2: Renormalization-group flows in the closed subspace 
H = Hb = for three cases of the inter-community cou- 
pling strength K/ J: (a) K/ J — 0.75, below the threshold 
value (K/J) m = ln(i(43 + 24V3))/ln(2 + V$) w 1.549; 
(b) K/J = {K/J) m ; (c) K/J = 3 > (K/J) m . In each 
case, the diagonal straight line represents the initial condi- 
tion J = Jq — 1/T, the thin vertical lines with arrows are 
sample flows, and the thick gray lines represent fixed points, 
with dark gray corresponding to stable fixed points, and light 
gray corresponding to unstable fixed points. The point where 
the light and gray curves meet is marginally stable. The 
dashed line marks the critical J c separating the disordered 
phase, flowing to the phase sink J* = 0, and the algebraically 
ordered phase, flowing to a line of finite temperature fixed 
points J* (Jo) (the dark gray curve). 

pairs satisfying 1^ = i. While G(£) cannot be di- 
rectly calculated from the RG flows, we can determine 
its long-distance scaling properties. Moreover, the aver- 
age correlation for a certain subset of pairs in the lattice 
can be explicitly calculated: at the nth level of the hi- 
erarchy, let site i be a hub of a community, and j be a 
site at the very edge of the same community, separated 
by a distance £ n = 2 11+1 — 1. Denote the average of Gij 
restricted to this subset as Ghub(^n)- After n RG steps, 
such pairs become nearest-neighbors along a solid 

bond, and thus we can obtain their thermodynamic av- 
erage through the conjugate recursion relation 36]. The 
number of such pairs is 2 2t ~ 2n ~ 1 . For 1 <C l n "C D, 
compared to the overall number of pairs with the same 
separation, Ng n ~ 2 3t , the subset forms a vanishingly 
small fraction of the total in the t — > oo limit. 



B. Phase diagram and critical properties 

Fig. [5] depicts various cases for the renormalization- 
group flows, and the corresponding phase diagram in 
terms of temperature T = 1 / J versus K/J at H = 
is shown in Fig. OHa). The two phases, both with M = 0, 
are (1 ) a disordered phase where pair correlations decay 
exponentially, G{1) ~ cxp(— £/£) with a finite correla- 
tion length £; (2) a phase with algebraic order, £ = oo, 
characterized by power-law decay of correlations, G(£) ~ 
£-v(t,k/J) _ Since this latter phase flows under RG to 
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FIG. 3: (a) Phase diagram. The black and gray curves repre- 
sent infinite- and second-order transitions respectively. (b,c) 
Critical exponents. The dotted line marks (K/J) m . 
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a line of finite-temperature fixed points (the dark gray 
curves in Fig. [3]), a different fixed point for every value 
of T and K/J, we have a varying exponent r/(T, Kj J). 
Fig. [3](b) plots 77 for K/J = 0.1 and 1, and we note 
that 77 — > for T — ► 0, as the system asymptotically 
approaches true long range order at T — 0. Strengthen- 
ing inter-community coupling has a similar effect, with 77 
decreasing for larger K/J. 

For K/J below a threshold value (K/J) m = 
In (yj-(43 + 24^3)) / In (2 + V5) « 1 -549, the phase tran- 
sition is infinite order: as T — > T c + we have exponential 
singularities of the Berezinskii-Kosterlitz-Thouless [HI, 
IH (BKT) form, just like in the XY model, f - e^/^ 
and the singular part of the specific heat C s j ng ~ e~ B ^, 
where t = (T — T c )/T c and the constants A,B > 0. It 
is interesting to note that for K/ J = 1 , our Hamilto- 
nian can be mapped by duality transformation to the 
Ising model on the small-world hierarchical lattice of 
Ref. [32| (the p = 1 lattice). On this dual network a 
similar BKT transition occurs, though with the alge- 
braic order in the high temperature phase (much as the 
Villain version of the XY model is dual to a discrete 
Gaussian model describing roughcning,with the low and 
high-temperature properties reversed [40 j ) . BKT singu- 
larities have also been observed in an Ising and Potts 
system on an inhomogeneous growing network [4ll l42l | . 
For K/ J > (K/ J) m , on the other hand, the phase transi- 
tion is second-order: £ ~ t~ u , C s i ng ~ t~ a for exponents 
a, v, plotted as a function of K/J in Fig.[3]Jc). Thus with 
increasing K/J the system almost looks like an ordinary 
second-order Ising transition: a critical T c below which 
we have something very close to long-range order, since 
77 is nearly zero. 



C. Griffiths singularities 

The disordered phase in our network differs in one im- 
portant aspect from a conventional paramagnetic phase: 



FIG. 4: (a) Magnetic susceptibility \ at K/J = 1 for T = 
2, 2.5 > T c . (b) x at K/J = 1 for T = 0.75, 0.25 < T c . (c,d) 
Hub correlation length £hub and specific heat C at K/J — 
1, 2.5. (d) Hub correlation function Ghub(^n), where £ n = 
2 n+1 - 1, at K/J = 1 for T = 0.90, 0.95, 0.98 < T c . 



M ~ H(l — In if) for small H, leading to a divergence in 
the susceptibility, \ ~ — \nH. We see this in Fig. [D^a) 
for T = 2, 2.5, and K/J = 1. As mentioned above, the 
mechanism for this nonanalyticity at H = is similar to 
the one behind the Griffiths singularity in random ferro- 
magnets. We can understand it through the distribution 
of Yang-Lee zeros [43[ of the partition function in the 
complex magnetic field plane. Introducing the variable 
z = e~ 2H , the Yang-Lee theorem states that the zeros of 
Z lie on the unit circle z = e in the complex z plane. If 
g(6) is the density of these zeros (a continuous distribu- 
tion in the thermodynamic limit), then for a regular fer- 
romagnet above T c we have g{9) =0 for a finite range of 
9 near the real axis 9 = , so that Z is analytic at H = . 
However in our case g{9) "pinches" the real axis: g(0) = 
but g(9 7^ 0) > 0. Since g{9) is related to the magnetiza- 
tion through g(9) = ^ lim r _ tl - ReM(z = re z9 ), we can 
deduce from the observed singularity that g(9) ~ 9 for 
small 9. The dominant contributions to this g(9) come 
from large communities, centered at hubs with degree 
k = 2 m + 1 for m ^> 1, despite their small probability 
P{k) = 3 • 4r m . 

We can see these contributions explicitly for the K/ J = 
case, adapting arguments used to derive scaling forms 
for g(9) in disordered ferromagnets [3 E^]. The system 
in this case is a disjoint set of solid-bond clusters, with 
the probability of a randomly chosen site being part of 
a cluster of size N m — 2 m ~ 1 +1 given by P m = 3(2 + 
2 m )/2 2m+1 (for m > 2). The average magnetization per 
site of such of cluster is easily calculated analytically, and 
takes the following approximate form for small H , 



M m « tanh(i/7 m (J)) , 



(4) 
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where 

fm(J) 

= 1 + N ™~ 1 tanh(2 J) (2 + (N m - 2) tanh(2 J)) . ^ 

The Hf m (J) term in Eq. (|4]) can be interpreted as 
the effective field felt by the cluster, with the function 
fm(J) varying between 1 at J = and N m at J = oo. 
In the large m limit f m (J) — > iV m tanh 2 (2 J) for all 
J > 0. To find g m (0), the cluster's contribution to 
the overall <?(#), we plug in a small complex magnetic 



field H = i(e 



corresponding to z = (1 — e)e l 



and take the real part of the resulting magnetization: 
g m (0) = ^ lim^o ReM m (H = ±(e - This gives 



5m (#) = lim 



sinh(e/ m (J)) 



2tt cos(0/ m (J)) + cosh(e/ m (J)) 



(G) 



This expression is dominated by high, narrow peaks at 
(9 = (2n + l)n/ f m (J), n = 0,1,..., and we can write it 
as a sum of delta functions in the small e limit, 



g m (0) « lim- V 

e— >0 7T ' — • 



efm(J) 



^ /m(J) 5 



n=0 



(2n + l)ir 

/m(J) 

Thus the total g{0) for the system is 



m— 2 n— 



(2n + l)7r 



(7) 



(8) 



For small 0, the nonzero contributions to <?(#) come from 
to values where f m (J) = (2n + l)7r/0 for some n. Since 
fm(J) ~ A r m tanh 2 (2J) for large to, these are the con- 
tributions of clusters with large size N m cx 1 /9, with 
a corresponding probability P m « 3/4iV m . Eq. (jHJ) be- 
comes 



v 30 2 tenhf(2J) / 



(2n + l)7T \ 



Nm tanh" 



(2 J) J 



(9) 



As to — > oo the delta function peaks become densely 
spaced, and from Eq. ([9]) it is evident for small 8 that 
g(6) scales like g(bd) = bg(9) for any constant b, con- 
sistent with the observation of g(ff) ~ 6* deduced from 
the singularity in M. Thus we see this behavior is 
directly related to the presence of large communities 
around highly connected hubs, which have a scale-free 
distribution P m ~ N^ 1 - 

In comparison, for bond-diluted ferromagnets below 
the percolation threshold large connected clusters are ex- 
ponentially rare, the resulting g{ff) ~ e~ A ^ T ^^, and the 
Griffiths singularity is much weaker, leading to a finite 
X at H — [III. Turning now to the algebraic phase 



for T < T c , here M and % behave as if at a critical 
point: M ~ ijV«(T,x/J) ) x „ h i/5(t,k/J)-i as # ^ . 

Fig. Hfb) shows x for IT/ J = 1, T = 0.75, 0.25, and 
Fig. |3I>) plots the exponent 5{T,K/J) for if/ J = 0.1 
and 1 . The corresponding scaling of the density of zeros 
isg(0) ~ 0V'Cr.*/J) n e a r = 0. 



D. Pair correlations 

Finally, we consider the behavior of the hub correla- 
tion function Ghub(^n)- For T > T c , we find an exponen- 
tial decay with l n , which we characterize by a correla- 
tion length £hub- The divergence of £hub as T — > T+ is 
shown in Fig. |3Jc) for K/ J = 1,2.5. Like the overall pah- 
correlation length £, ^hub diverges with a BKT form for 
K/J < (K/J) m and as a power law for K/J > (K/J) m . 
The onset of the rapid increase in £hub coincides with the 
position of the peak in the specific heat C, plotted in 
Fig. QJd) . Just like in the XY model [33], C is smooth 
at T c for all K/J, and the peak occurs at T > T c , corre- 
sponding to the onset of short-range order in the system. 
For T < T c , Ghub(^n) has a surprising behavior: as seen 
in Fig.0|e), it approaches a nonzero limit 
signature of long-range order. However, since Ghub(^n) 
describes only a subset of pairs, a vanishingly small frac- 
tion of the total for large £ n , the long-range ordering of 
these pairs is compatible with M being zero. The pres- 
ence of such long correlations in the algebraic phase, and 
the overall slow power-law decay of G(£) with I, is re- 
markable given that for "fat-tailed" scale-free networks 
(i.e. with 7 < 3) pair correlations longer than nearest- 
neighbors are typically suppressed: one can prove that 
G(£ > 1) = at H = in the thermodynamic limit, 
if x(H = 0) is finite [3^]. In our case x(if = 0) = oo 
at all T, the proof does not apply, and we see that this 
fractal modular lattice is an important exception to the 
general expectation of weakened pair correlations on net- 
works [al. 



IV. CONCLUSIONS 

In conclusion, we have introduced a hierarchical lat- 
tice network with the modular structure and fractal scal- 
ing characteristic of a wide array of real-world networks. 
The Ising model on this lattice — solved through an ex- 
act RG transformation — exhibits an interesting transi- 
tion. A disordered phase with Griffiths singularities gives 
way at low temperatures to algebraic order: the system 
behaves as if at criticality for a broad range of param- 
eters, and we find power-law decay of pair correlations, 
unexpected for this type of scale-free network. 

The thermodynamic phenomena observed here are not 
confined to one particular network. In fact, we can con- 
sider the much larger class of hierarchical lattices that 
form scale-free networks on which the Ising model ex- 
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hibits a standard order-disorder transition: these include 
fractal lattices on which Migdal-Kadanoff recursion rela- 
tions are exact [H, , related hybrid lattices [Ho[ , and 
their duals [HJ . In all these cases we can modify the con- 
nected graph which defines the lattice construction step 
as follows: replace a subset of the bonds in the graph 
with dashed bonds (which remain unaltered as we iter- 
ate the construction) in such a way that the graph would 
break into two or more disjoint pieces if the dashed bonds 
were cut. The Ising model on the resulting hierarchical 
lattice will no longer flow under renormalization to an 
ordered fixed point at low temperatures, but rather to a 
continuous line of fixed points, yielding an algebraically 
ordered phase. And the power-law distribution of highly 
connected hubs in such networks will lead to Griffiths sin- 
gularities in the disordered phase. Conversely the duals 



of such networks — which have infinite fractal dimension 
and show small-world scaling — have algebraic order at 
high temperatures. Recent studies have highlighted the 
diversity of structural properties in families of hierarchi- 
cal lattice networks [3J, y4j — the ability to tune degree 
exponents, fractal dimensionality, and other topological 
aspects of these networks by varying the defining graph. 
The structural richness of these networks is manifested 
through unusual phase transitions and critical phenom- 
ena, already apparent even in a simple system like the 
Ising model. The use of renormalization-group methods 
to characterize cooperative behavior on this broad class 
of networks, both in the pure case and in the presence of 
bond randomness [32j] , will be the subject of future work. 

I thank A.N. Berkcr, T. Garel, and H. Orland for useful 
discussions. 
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